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ABSTRACT 

Context. Broadband secondary-eclipse measurements of strongly irradiated hot Jupiters have indicated the existence of atmospheric 
thermal inversions, but their presence is difficult to determine from broadband measurements because of degeneracies between molec¬ 
ular abundances and temperature structure. Furthermore, the primary mechanisms that drive the inversion layers in hot-Jupiter atmo¬ 
spheres are unknown. This question cannot be answered without reliable identification of thermal inversions. 

Aims. We apply high-resolution (R = 100000) infrared spectroscopy to probe the temperature-pressure profile of HD 209458 b. This 
bright, transiting hot-Jupiter has long been considered the gold standard for a hot Jupiter with an inversion layer, but this has been 
challenged in recent publications. 

Methods. We observed the thermal day side emission of HD 209458 b with the CRyogenic Infra-Red Echelle Spectrograph (CRIRES) 
on the Very Large Telescope during three nights, targeting the carbon monoxide band at 2.3 |im. Thermal inversions give rise to 
emission features, which means that detecting emission lines in the planetary spectrum, as opposed to absorption lines, would be 
direct evidence of a region in which the temperature increases with altitude. 

Results. We do not detect any significant absorption or emission of CO in the day side spectrum of HD 209458 b, although cross¬ 
correlation with template spectra either with CO absorption lines or with weak emission at the core of the lines show a low-significance 
correlation signal with a signal-to-noise ratio of ~ 3 - 3.5. Models with strong CO emission lines show a weak anti-correlation with 
similar or lower significance levels. Furthermore, we found no evidence of absorption or emission from H 2 O at these wavelengths. 
Conclusions. The non-detection of CO in the day side spectrum of HD 209458 b is interesting in light of a previous CO detection in 
the transmission spectrum. That there is no signal indicates that HD 209458 b either has a nearly isothermal atmosphere or that the 
signal is heavily muted. Assuming a clear atmosphere, we can rule out a full-disc dayside inversion layer in the pressure range 1 bar 
to 1 mbar. 

Key words. Planets and satellites: atmospheres - Infrared: planetary systems - Methods: data analysis - Techniques: spectroscopic 
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1. Introduction 


1.1. Hot-Jupiter atmospheres 


Hot Jupiters form an ideal starting point for exoplanet atmo¬ 
spheric research. While they experience high levels of stellar 
irradiation, their atmospheric scale heights are relatively large, 
and their high temperatures result in high planet-to-star contrast 
ratios, especially in the infrared. Furthermore, since hot Jupiters 
have short orbital periods (1-5 days), they have a high probabil¬ 
ity of transiting their host star, and transits and eclipses occur on 
conveniently short timescales. The accessibility of hot Jupiters 
have made them excellent test objects for developing techniques 
for atmospheric characterisation - which are applicable to the 
next generation of telescopes and instruments and hence also 
to sm aller planets (e.g. [Schneider 1994[ Webb & Wormleatonj 


2001 [S^ger & Deming|20T0f Snellen et al.||2013| ). The com- 

position and structure of individual hot-Jupiter atmospheres are 
themselves of great interest to better understand the complex 
physical processes involved (including strong stellar irradiation. 


Based on observations collected at the ESO Very Large Telescope 
(Program 186.C-0289) 


atmospheric circulation, and photochemical processes), and be- 
cause they may hold clues to the planetary for mation history (e.g. 
Qberg et ^2011 Madhusudhan et al.|2014 ). 


Atmospheric characterisation of hot Jupiters has so far 
mainly focused on transiting planets because they offer unique 
observational possibilities. When a planet transits its host star, 
part of the stellar light is filtered through the planetary atmo¬ 
sphere, leaving an absorption fingerprint in the observed light 


([Seager & Sasselov|2000t Brown et al.|2QQl[[Charbonneau et al. 


|2002| ). At the opposite side of the orbit, the star occults the 
planet, providing an opportunity to isolate the planetary contri¬ 
bution to the total flux. The difference between observations of 
the combined light of the star and planet in and out of secondary 
eclipse can be used to determine the planet’s dayside spectrum 
( [Charbonneau et al.|2005[|Deming et al.|2005| ). In addition, flux 
variations during the entire orbit can be measured as a function 
of phase, ga ining information on the day- and nightsi de contri¬ 
butions (e.g. [Knutson et al.|20Q7at[Snellen et al.|2009| ). 


Hot Jupiters are assumed to have atmospheres that are dom¬ 
inated by molecular hydrogen, but the spectrum of a given 
planet will contain several spectroscopically active trace gases 
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of carbon- and oxygen-bearing molecules. At high temperatures 
and observed in the infrared, the most abundant of these mole¬ 
cules are expected to be CO and H 2 O, followed by CH 4 . Other 


significant molecules are possibly CO 2 , C 2 H 2 , and HCN (|Mad- 


|husudhan|2012|[Moses et al.|2013| ). 


1.2. Thermal inversion layers 

The upper atmospheres of hot Jupiters are thought to be radi¬ 
ation dominated, with the radiation efficiency depending criti¬ 
cally on the emission and absorption of the specific molecules 
that are present. Generally, the radiation from the host star is 
expected to deeply penetrate the atmosphere and heat it from be¬ 
low, which results in temperatures that decrease with increasing 
altitude ( |Guillot|2010| ). However, strong optical or UV absorbers 
in the upper layers of the atmosphere can cause thermal inver¬ 
sion - a layer or region in the atmosphere in which the tempera¬ 
ture instead increases with altitude ( jde Pater & Lissauer|2010| ). 
Thermal inversions are common among the solar system plan¬ 
ets. Jupiter, along with the other giant planets, has thermal inver¬ 
sions caused by CH 4 -induced hazes, while the inversion in the 
terrestrial stratosphere is caused by O 3 (e.g. |Seager||2010| and 
references therein). 

Thermal inversions may be commonly present in hot-Jupiter 
atmospheres as well. Secondary-eclipse measurements with the 
Spitzer Space Telescope have made it possible to probe the ver¬ 
tical temperature structure of bright transiting hot Jupiters. By 
measuring the relative depths of a secondary eclipse in multi¬ 
ple band-passes, a low-resolution thermal spectrum of the exo¬ 
planet can be constructed that can then be compared with the¬ 
oretical spectra. There have been several reports of inversion 
layers in hot Jupiters from Spitzer observat ions (HD 209458b, 
Burrows et al.||2QQ7l [Knutson et al.||2008t XO- lb, [Machalek 


et al.|2008t XQ ^-2b, Machalek et al.|2069t TrES-4, [Knutson et al. 
2 OO 9 I TrES-2, [Q’Donovan et al.[[2010| Hat-P-7b, [Christiansen 


(e.g. HD 189733b, [Charbonneau et al.|2008[ ). The inference of a 
thermal inversion is in most cases based upon an excess of flux 
in the 4.5 [im and 5.8 [im Spitzer band-passes, which has been 
interpreted as due to water emission features. However, many 
of these claims are based upon the comparison of only a few 
inverted and non-inverted atmospheric models, which may not 
adequately map degeneracies between atmospheric parameters. 
In contrast, a systematic retr ieval analysis of s econdary-eclipse 
spectra of nine hot Jupiters ( [Line et aT][2Q14| ) found little evi¬ 
dence for thermal inversions over a wide range of effective tem¬ 
peratures (with the exception of HD209458b). 

The nature of the possible responsible absorbers is as yet un¬ 
known in this high-temperature regime. It has been proposed 
that highly irradiated planets are warm enough for significant 
amounts of gaseous TiO and VO to exist in the upper atmo¬ 
sphere, providing the necessary opacity to generate a tempera- 
ture inversion ([Hubeny et al.[[2003[ [Fortney et al.[[2006] [2008[ 
[Burrows et al.|2008[ ). This advocates a correlation between high 
irradiation and existence of inversion layers. No observational 
evidence for TiO and VO absorption in hot Jupiter atmospheres 
has been found, however ([Huitson et al.|2013t[Sing et al.|2013 


[Gibson et al.|2013|[Bento et al.|2014t[Hoeijmakers et al.|2Q14[ ), 

although [Desert et al.[p008[ ) suggested TiO and VO molecules 
as possible candidates to explain observed absorption broadband 
features in the optical transmission spectrum of HD 209458 b. 

The TiO/VO hypothesis has since been challenged by 
[Spiegel et al.[ ( [2009] ), who argued that even a ten times supersolar 
abundance of VO is insufficient to drive a thermal inversion and 
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that the high molecular weight of TiO requires substantial verti¬ 
cal mixing to retain a high abundance of TiO at low pressures. 
Apart from a deep vertical cold trap, TiO can be depleted by a 
nightside cold trap if it condenses to particles la rger than a few 
micro ns ([Parmentier et al.[[2013| ). Furthermore, [Madhusudhan[ 
[et al.[ ( [20il| ) found that for carbon-rich (C/0 >1) hot Jupiters, 
most of the oxygen is bound up in CO, and as a result, the abun¬ 
dances of TiO and VO are too low to cause a thermal inversion. 
Alternatively, [Knutson et al.[ ( [2010| ) found a preference for plan¬ 
ets thought to host an inversion to orbit around chromospheri¬ 
cally quiet stars and vice versa, and they theorised that increased 
UV levels for active stars were destroying the unknown com¬ 
pounds responsible for inversion layers. The inversion-causing 
absorbers could also be the result of non-equilibrium processes. 
[Zahnle et al.[ ( [2009[ ) propounded the photochemically produced 
sulfur compounds, S 2 and HS, as efficient UV and optical ab¬ 
sorbers at high temperatures. 

Comparative studies of hot Jupiters with and without temper¬ 
ature inversions are key to explaining the physical causes of in- 
version layers. Although such studies have been attempte d ( |Fort-[ 


et al.[[ 2010 | , while other planets show no sign of an inversion 


ney et al.|2008[ [Burrows et al.|2008[ [Knutson et al.|2Q10[ ), this is 
not truly possible until a reliable method is available for measur¬ 
ing thermal inversions. Broadband photometry does not resolve 
the molecular spectral bands, which means that emission lines 
cannot be directly detected, and the inference of a thermal in¬ 
version must therefore rely heavily on models and their inherent 
assumptions. The Spitzer inferences of inversion layers tend to 
rely on assumptions about the presence of water, and typically, 
solar composition is assumed. [Madhusudhan & Seager[ ( [2010[ ) 
explored a wider paramater space and found evidence for degen¬ 
eracies between the temperature pressure profile and the molec¬ 
ular abundances. 


1.3. High-resolution spectroscopy 

Ground-based high-dispersion spectroscopy in the near-infrared 
has recently become successful in characterising hot Jupiters. By 
resolving molecular features into individual lines, the lines can 
be reliably identified through line matc hing ([Brown et al.|200^ 
[Deming et al.[2005] Barnes et al.[2007[ ). [Snellen et al.[ ( |201Q[ ) em- 
ployed the Cryogenic Infra-Red Echelle Spectrograph (CRIRES, 
[Kaufl et al.|2004[ ) at the Very Large Telescope (VLT) to observe 
the transmission spectra of HD 209458 b with a resolution of 
R = 100000, and detected CO at 2.3 qm. The technique makes 
use of the changing radial velocity of the planet during the ob¬ 
servations to separate the planetary lines from the telluric and 
stellar lines. 

The target planet can be observed favourably either during a 
transit (transmission spectroscopy) or at phases shortly before or 
after superior conjunction (dayside spectroscopy). Since a sec¬ 
ondary eclipse is not required to separate the contamination from 
the parent star, dayside spectroscopy can also be applied to non¬ 
transiting planets, r Bootis b has been observed at 2.3 [im with 
high resolution by [Brogi et al.[ ( 20 T^ and [Rodler et al.[p 012 | ), 
who detected absorption of CO, and the radial velocity mea¬ 
surements led to the first mass determination of a non-transiting 
hot Jupiter. The dayside spectrum of 51 Peg b revealed absorp¬ 
tion from CO and possibly water ( [Brogi et~ar][2013[ ), but the 
signal was not detected in one of the three nights of observa¬ 
tions, and therefore more data are necessary to confirm the re¬ 
sult. CO absorpti on was detected in 2. 3 qm dayside observations 
of HD 189733 b ([de Kok et al.|2013[). H 2 O absorption was de¬ 
tected in 3.2 [im dayside observations of the same planet ( [Birkby 
[et al.[[2013] ), demonstrating the ability of high-resolution spec- 
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troscopy to also detect a more complex molecule in a wavelength 
regime more heavily contaminated by telluric lines. Recently, a 
combined signal from CO and H 2 O absorption was detected in 
2.3 [im dayside observations of HD 179949 b ( [Brogi et al.|2014| ), 
and|Lockwood et al.|(|2014|) detected H 2 O in the dayside spec- 
truni of T Bootis b at L-band using NIRSPEC at Keck. 

Dayside high-dispersion infrared spectroscopy has the po¬ 
tential to provide evidence for an inversion layer in a hot Jupiter. 
The molecular bands are resolved into individual lines, so that 
the average shape of the lines and the average contrast with re¬ 
spect to the stellar continuum can be measured, or at least con¬ 
strained. The dayside spectrum of a hot Jupiter is sensitive to 
the presence of thermal inversions because the line-of-sight to 
a high degree coincides with the vertical direction in the plan¬ 
etary atmosphere. The sign of the vertical temperature gradient 
determines whether the thermal emission from the molecules in 
the upper atmosphere is seen as emission or absorption relative 
to the continuum part of the spectrum, emitted deep in the at¬ 
mosphere (between 1 and 0.1 bar). The sign of the measured 
line contrast thus allows distinguishing between emission and 
absorption lines, and the shape of the lines can help further con¬ 
strain the temperature-pressure prohle. 


1.4. HD 209458 b 

At this point, ah published high-dispersion detections of mole¬ 
cules have detected lines in absorption, but if instead emission 
lines were to be detected, this would prove the existence of a 
thermal inversion. We therefore apply here the high-resolution 
technique to the dayside of the hot Jupiter HD 209458 b. The 
planet was the hrst exoplanet known to transit its parent star 
( Charbonneau et al.|2000| ), and today it is one of the best-studied 
hot Jupiters. 

In the context of inversion layers, HD 209458 b is of par¬ 
ticular interest since it was the hrst exoplanet reported to have a 
thermal inversion ( [Burrows et al.|2007l[Knutson et al.|2008] ), and 
until recently, it was often proclaimed to be the benchmark for 
a hot Jupiter with a thermal inversion in its upper atmosphere. 
This was based on secondary-eclipse observations at 3.6, 4.5, 
5.8 and 8.0 [im usin g the Spitzer Infrared Array Camera (IRAC 
Knutson et al. 12008] ). The 4.5 and 5.8 qm bandpasses showed an 


excess of hux relative to the 3.6 qm bandpass, in contrast to tra¬ 
ditional models that expected a trough at these wavelengths due 
to water absorption features. The excess hux was best explained 
with a low-pressure inversion layer producing water emission 
features. The claim of a thermal inversion in HD 209458 b was 
supported by extensive modelling by [Madhusudhan & Seager[ 
( [2009[[20T0| and [Line et al.[ ( [2014| ). However, the existence of a 
thermal inversion in the atmosphere of HD 209458 b has recently 
been challenged. [Zehem et al.[ ( [2014[ ) measured the full orbit 

4.5 [im phase curve with Spitzer IRAC using the now standard 
staring-mode observations and pixel-mapping techniques. They 
found the 4.5 qm secondary-eclipse depth to be 35% shallower 
than previously measured by [Knutson et al.[ ( [2008[ ), who did not 
use staring mode, nor intra-pixel sensitivity maps. This revision 
does not rule out an inversion layer, but is consistent with models 
both with and without a thermal inversion, or with a blackbody. 
[Diamond-Lowe et~aL] ( [2014[ ) have performed a self-consistent 
analysis of all available Spitzer/IRAC secondary-eclipse data of 
HD 209458 b. This includes both a reanalysis of the eclipses in 
the four IRAC bandpasses published by [Knutson et al.[ ( [2008[ ) 
and hve previously unpublished eclipses in staring mode, one 
in 2007 at 8.0 qm, two in 2010 at 4.5 [im, and two in 2011 at 

3.6 [im. They found no evidence for a thermal inversion in the 


atmosphere of HD 209458 b. Their best eclipse depths are well 
htted by a model where the temperature decreases between pres¬ 
sure levels of 1 and 0.01 bars. 


1.5. Re-evaluation of the HD209458b CO abundance from 
Snellen et a/.[{2oiop 


Carbon monoxide absorption has previously been observed with 
CRIRES at the VET during a transit of HD 209458 b ( [Snehen[ 
et al.[[MT0| , with a relative line strength at the level of 1 - 
1.5 X 10“^ with respect to the host star spectrum. Although a 
constraint on the CO volume-mixing ratio was presented, in 
hindsight it became clear that the underlying analysis did not 
include a thorough investigation of the parameter space. Eur- 
thermore, an erroneous conversion factor in the computation of 
the pressure broadening of the absorption lines resulted in a too 
high derived abundance. We note that the conversion error was 
only present in [Snellen et al.[ ( [20T0| ) and not in any of our subse¬ 
quent work. We since performed a re-analysis of the constraints 
to the CO volume mixing ratio that can be inferred from the 
high-dispersion transit observations. A large grid of models was 
constructed, assuming a range of temperature-pressure prohles, 
abundances of a range of molecules, and pressure levels of the 
continuum. This re-analysis indicates that the observed relative 
line strength corresponds to a CO volume-mixing ratio in the 
range 10“^ - 10“'^, with an additional uncertainty of a factor of 
~2 if uncertainties in the high-altitude atmospheric temperature 
at the terminator region are included. However, if clouds or hazes 
contribute to the continuum extinction, the corresponding CO 
abundance could in turn be signihcantly higher. 


1.6. Outline 

In this paper we present three nights of CRIRES observations 
of HD 209458 at 2.3 qm, targeting CO in the planet’s dayside, 
allowing us to probe the vertical temperature structure of the at¬ 
mosphere. Section 1^ details our observations of HD 209458 b, 
and Sect, [^describes the data analysis steps from the raw spec¬ 
tra to removal of telluric contamination. In Sect. [4] the cross¬ 
correlation analysis we used to extract the potential planetary 
signal is explained, and the results are presented in Sect. We 
discuss these in Sect.[^ and conclude in Sect.[^ 


Table 1. Details of the observations, showing the observing date for the 
beginning of the night in local time, the use of nodding, the planetary 
phase, the total observing time, the number of observed spectra, range 
in airmass, water vapour, detector integration time (DIT), number of 
integrations per nod or observation (NDIT), and detector counts. 



2011-08-04 

2011-09-05 

2011-09-12 

Nodding 

yes 

yes 

no 

Orbital Phase 

0.51-0.57 

0.55-0.62 

0.54-0.61 

lobs [h] 

6.07 

5.53 

6.02 

Nobs 

118 

108 

122 

Airmass 

2.16-1.38 

2.27-1.38 

2.20-1.38 

DIT [s] 

30 

30 

50 

NDIT 

5 

5 

3 

Cnts [#/pix/s] 

40-65 

50-70 

50-90 
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2. Observations 


We observed HD 209458 (K = 6.31 mag) for a total of 17.5 
hours during three nights in August and September 2011 as part 
of the large ESO program 186.C-0289, which was designed to 
study the brightest transiting and non-transiting hot-Jupiter sys¬ 
tems visible from Cerro Paranal. The observations were obtained 
with the CRIRES instrument, which is mounted at the Nasmyth 
A focus of VET Antu. We observed at the standard wavelength 
settings for the reference wavelength 2.3252 [im, covering the 
ro-vibrational (2,0) R-branch of carbon monoxide. The CRIRES 
instrument has four Aladdin III InSb detectors, each of the size 
1024x512 pixels with gaps of about 280 pixels between the in¬ 
dividual detectors. 


We used a 0.2" slit in combination with the Multi Ap¬ 
plication Curvature Adapative Optics system (MACAO, 


Arse- 


|nault et JT]|2003| ) to achieve the highest possible resolution of 
R ^ 100000. During each of the three nights, the target was ob¬ 
served continuously for 5-6 hours at planetary phases shortly af¬ 
ter a secondary eclipse, when a significant portion of the dayside 
of the planet was visible. Additionally, a standard set of calibra¬ 
tion frames was obtained during daytime and twilight. 

The first two nights were observed in nodding mode, where 
the telescope is nodded along the slit by 10" in an ABBA pat¬ 
tern, providing an easy method for background subtraction. The 
third night was observed without nodding in an effort to enhance 
stability and reduce overheads. Details of the observations are 


given m 


99 


3. Data analysis 


The data in this work have been analysed in a similar 
in [Snellen et al.| ( [2010| ) and jBrogi et al.j ( [20121 2013[ ), 
analyses were performed with new purpose-built Python 


way as 
but all 
scripts. 


3.1. Extracting the one-dimensional spectra 

The basic image processing was performed with the CRIRES 
pipeline version 2.1.3 and the corresponding version 3.9.0 of 
ESOREX. The images were dark-subtracted and fiatfielded and 
were corrected for bad pixels and non-linearity effects. Eor the 
first two nights, the CRIRES pipeline was also used to combine 
the images in AB nodding pairs, thus performing a background 
subtraction, and to extract the ID spectra. This resulted in a sin¬ 
gle extracted spectrum for every two observations. These spec¬ 
tra were extracted with the optimal extraction technique ( [Home 
[1986[ ). The spectra of the final night were observed without nod¬ 
ding, for which the CRIRES pipeline was ill-suited. Instead, we 
used the IRAE data analysis package apaii, both for the back¬ 
ground subtraction and the optimal extraction of the spectra. 

Detectors 1 and 4 are both affected by the odd-even effect, a 
non-linear change in gain between odd and even columns along 
the spectral directioij^ This is the result of the reading direc¬ 
tion being perpendicular to the spectral dispersion for these two 
detectors. In particular, detector 4 showed residual odd-even im¬ 
prints after non-linearity correction. Eurthermore, detector 1 has 
very few expected carbon monoxide lines. Eor these reasons, we 
chose to leave out detectors 1 and 4 from further analysis, which 
now covers the wavelength range 2.3038 [im-2.3311 qm. 


step 1 


step 2 


■ step 3 



. - : - • . ,. 

step 4 

—<— 1 —j—t 1. j—|. r ■ii“r ..— ] — i —t—h~ 3—1 ■ 1 - i—— i —J -I -1 - j 

■ step 5 

1 1 j r r L J |. 1 ] -J. .1 1 X \ 1 -t ] . L. J .i.. I J J .1 r. 1 ..j 


2322 2322.5 2323 2323.5 2324 2324.5 


1.04 

1.03 

1.02 

1.01 

1.00 

0.99 

0.98 

0.97 

0.96 


Wavelength [nm] 


Fig. 1. Overview of the different steps in the data analysis, showing 
a small fraction of the third array from the first night of observations. 
From top to bottom, step 1 shows the extracted spectra after bad-pixel 
correction, step 2 after spectral alignment and normalization, step 3 after 
first-order airmass correction, step 4 after second-order airmass correc¬ 
tion, and step 5 after normalization of each column with the variance. 


3.2. Bad-pixel correction and wavelength calibration 


The extracted spectra were treated separately per night and per 
detector. The spectra were arranged into two-dimensional ma¬ 
trices with wavelength along the x-axis and time along the y- 
axis. The bad-pixel correction included in the CRIRES pipeline 
is insufficient, and we therefore manually identified bad pixels, 
columns, and regions through visual inspection of the matrices 
with the program DS9. A few larger regions consisting of several 
adjacent bad columns, such as the columns closest to the edge of 
the CCD’s and a well-known defect on the second detector, were 
masked during the subsequent data analysis. The masked pixels 
made up about 2% of the total number of pixels. Remaining in¬ 
dividual bad pixels and bad columns were corrected with cubic 
spline interpolation based on the four nearest neighbours on each 
side in the row. 


All spectra of a given matrix were aligned to a common 
wavelength grid. We determined the offsets between centroids 
(determined using a Gaussian fit) of deep and isolated telluric 
lines in the individual spectra (typically 15-20 per array) and 
those in the average spectrum. The centroid offsets from a given 
spectrum showed no trend with wavelength and a typical scat¬ 
ter of less than 0.1 pixel with respect to each other. Each spec¬ 
trum was therefore shifted with a spline interpolation using the 
median centroid offset of that spectrum. Next, we matched the 
pixel positions of centroids of telluric lines in the average aligned 
spectrum with the wavelengths in a synthe tic telluric transmis¬ 
sion spectrum from ATRAhj^ (Lord 1992). The pairs of pixel 
position and wavelength were fitted with a second-order poly¬ 
nomial to yield the wavelength solution for a given night and 
detector. The highest residuals to the second-order fits were of 
the order +3 x 10“^ qm, which corresponds to 20% of a pixel. 
Einally, each spectrum was normalised by its median value. 


^ http://www.eso.org/sci/facilities/paranal/instruments/crires/docWLT- _ 

MAN-ESO-14500-3486_v93.pdf ^ http://atran.sofia.usra.edu/cgi-bin/atran/atran.cgi 
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3.3. Removing telluric contamination 

The observed near-infrared spectra are dominated by absorp¬ 
tion lines originating in the atmosphere of Earth, and a crucial 
datareduction step is therefore to remove this telluric contami¬ 
nation. The telluric lines are stationary in wavelength (but vary 
in strength) over the course of the observations, and therefore fall 
on fixed columns in the matrices. On the other hand, molecular 
lines arising in the planetary atmosphere are Doppler-shifted in 
each consecutive spectrum as the radial component of the or¬ 
bital velocity of the planet changes, resulting in diagonal fea¬ 
tures across the matrices, buried in the noise. The planet lines 
typically move by 30 to 40 pixels during one night of observa¬ 
tion. This means that we can manipulate the matrices column by 
column without significantly affecting the planetary signature. 

The variations along columns of the matrices are strongly 
correlated with the changing airmass over the time interval of the 
observations. We corrected for this with a linear least-squares fit 
of the logarithm of the fluxes in each column as a function of 
geometric airmass. 

At the positions of some of the more prominent telluric lines, 
there are clear residuals, showing a correlation with time. These 
second-order effects stem from changing conditions in Earth’s 
atmosphere over the course of the night. Two distinct patterns in 
time were recognizable and were found to correspond to water 
and methane lines. We sampled the two effects of flux with time 
in a few affected columns in a strong water line and a strong 
methane line. We then corrected the remaining matrix by ex¬ 
pressing each column as a linear function of the sampled effects 
and dividing by the fit. The second-order airmass fitting can be 
affected by the planetary signal we are ultimately looking for. In 
particular, in columns (wavelengths) where the planet signal is 
at the start or end of the time series, the fit can be affected by 
this additional signal and therefore partially compensate for it. 
This is overall typically a 10% effect. Therefore we masked pix¬ 
els expected to contain a planetary signal when performing this 
step. 

Eollowing this, a few columns that coincided with the wave¬ 
lengths of the deepest telluric lines were visibly more noisy than 
the rest. We divided each column by the variance of itself to 
ensure that the noisiest columns did not dominate the cross¬ 
correlation in the next step of the data analysis. An example of 
the intermediate results of the different steps in the data reduc¬ 
tion is shown in Eig.[2 We achieve a typical signal-to-noise ratio 
(S/N) of -250 per pixel per spectrum for the first two nights, and 
-200 per pixel per spectrum for the final night, for which we did 
not dither. 

4. Search for the planet signal 

The individual planetary lines are buried in the noise of the resid¬ 
ual spectra. With the chosen wavelength setting of 2.3 qm, we 
are targeting the (2,0) R-branch of CO, and we expect tens of 
CO lines within the observed wavelength range (see Eig.|^. The 
signals from the individual lines are combined by cross correlat¬ 
ing the residual spectra one at a time with purpose-built model 
spectra. In this section we describe the model spectra and the 
cross-correlation analysis used to search for a planetary signal in 
the observed spectra. 

4.1. Model spectra 

The model spectra of HD 209458 b were calculated line by line 
using parameterised temperature-pressure profiles (T/P profiles). 


H 2 -H 2 collision-induced absorption (|Borysow et al.||200I[ |Ba 
|rysow||2002I) and with either CO or Ho O as a single trace gas 


gas. 


rH 20 data were taken from HITEMP 2010 ([Rothman 


et al.|2010| ), and a Voigt line profile was employed. We used the 


parameter values for the Voigt profile as given by the HITEMP 
database. Uncertainties in the line shape parameters will mainly 
influence the depths of the lines, not the cross-correlation values, 
since the CO lines are widely separated, and the spectral shape 
will not be affected much. In this article, we focus on the analy¬ 
sis and results from cross correlating with the model spectra with 
CO as the single trace gas. A short summary of the results from 
the H 2 O models, which showed no signs of a signal, is given in 
Sect.r 


We performed the cross-correlation analysis for a grid of 
models spanning a range of T/P profiles and CO volume-mixing 
ratios (VMR). The T/P profiles describe the average vertical tem¬ 
perature structure of the dayside of the planetary atmosphere. All 
the models in the grid are isothermal with a temperature of Ti = 
[1500 K, 2000 K] at pressures higher than Pi = 1 bar. The tem¬ 
perature decreases with a constant lapse rate (i.e. the rate of tem¬ 
perature change with log pressure) until it reaches T 2 = 1000 K 
at pressure P 2 = [10“^ bar, 10“^ bar]. At higher altitudes there 
is either an inversion layer or not, depending on the temperature 
T 2 = [750 K, 1750 K] at the pressure P 3 = [ 10“^ bar, 10"^ bar], 
above which the models again become isothermal. This amounts 
to eight different T/P profiles with an inversion layer and eight 
T/P profiles without, all of which were tested with three different 
volume-mixing ratios, VMR = [lO”"^’^, 10“"^'^, 10“^'^]. 

The three bottom rows of Pig. show representative exam¬ 
ples of CO model spectra from the T/P grid: A model atmosphere 
without a thermal inversion, and two with inversion layers, all 
with CO VMR = The model without an inversion has 

absorption lines, while the two models with thermal inversions 
conversely show emission at the core of the CO lines and absorp¬ 
tion in the wings. In the following, we refer to the model with 
the strongest emission lines as the strong-inversion model and to 
the one with weaker emission as the weak-inversion model. The 
distinction between weak and strong is based on whether or not 
the core emission is below or above the continuum level after 
convolving to the resolution of the CRIRES instrument. All the 
models from the grid without an inversion layer have relatively 
broad CO absorption lines (PWHM - 10 to 45kms“^), which 
have a tendency to smear out the planetary signal in the cross¬ 
correlation analysis. We therefore also investigated a CO model 
with more narrow absorption lines (PWHM -4 km sec“^), il¬ 
lustrated in the top row of Pig. This model has no thermal 
inversion and an adiabatic lapse rate until the assumed lowest 
temperature of 500 K, thus maximising the temperature gradi¬ 
ent for the given CO volume-mixing ratio, which in this case is 
1 X 10“^^ 


4.1.1. More inversion models 


Por comparison, we calculated an additional suit of CO mod 
els with a thermal inversion with T/P profiles qualitatively and 
quantitatively similar to the best-fit T/P profiles in Madhusud- 
han & Seager|(|2009'|). We call these the MS models. Madhusud- 
han & Seager| ( [200^ applied an atmospheric retrieval method 


with a parametric T/P profile to the Spitzer IRAC secondary- 
eclipse depths in [Knutson et al.j ( [2008[ ), supporting and con¬ 
straining a dayside thermal inversion in HD 209458 b. The T/P 
profiles of the MS models are shown in black in Pig. to¬ 
gether with the profiles from the thermal inversion models of 
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Fig. 2. Examples of CO model spectra used in the cross-correlation analysis. The model in the top row has an adiabatic lapse rate and a CO VMR 
of 10""^ ^, resulting in strong, narrow absorption lines. The following three rows are models from the T/P grid described in Sect. 4.1 all with CO 
VMR = 10“^ ^: the first is without a thermal inversion and with relatively broad absorption lines, the second is with a weak inversion, and the 
last model is for a strong thermal inversion. The left column displays the models before convolution to CRIRES resolution over the wavelength 
coverage of the two central detectors. The middle column zooms in on the CO line at approximately 2.3205 |im, showing the line profiles both 
before (solid) and after (dashed) convolution to CRIRES resolution. The right column gives the corresponding T/P profiles. 


the T/P grid described in Sect. |4.1| shown in red. We calculated 
the MS models using three different CO volume-mixing ratios, 
VMR = [10■^ 10“^, 10“^]. 


The thermal inversions of the MS models are located deeper 
in the atmosphere than those from the T/P grid in the previous 
section. This results in strong emission in all cases and in line 
shapes that are very sensitive to the volume mixing ratio. The 
models with the lowest VMR = 10“^ probe relatively low alti¬ 
tudes, making the inversion layer more spectroscopically domi¬ 
nant, thus resulting in line profiles with pure emission after con¬ 
volution to CRIRES resolution. Models with higher VMRs still 
have siginificant emission due to the deep inversion layer, but 
they have absorption at the core of the lines because the core 
probes higher altitudes than for the lower VMRs where the tem¬ 
perature again decreases with increasing altitude. The shape of 
the lines at CRIRES resolution of MS models with a VMR of 
10 “"^ or 10“^ resemble the line profile of the weak-inversion 
model in the middle panel of Eig.j^at CRIRES resolution, only 
mirrored vertically. 


Note that recent analyses of the dayside spectrum of 
HD209 458b by |Diamond-Lowe et al.| ( [MT4| ) and |Zellem et al.| 
( |2014| ) found that the data were best fit by a near-isothermal T/P 
profile. However, with high-dispersion spectroscopy we are only 
sensitive to the narrow components of the spectra, which are pro¬ 
duced relatively high up in the atmosphere, making our observa- 
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tions insensitive to such a near-isothermal atmosphere. This is 
further discussed in Sect, [b] 


4.2. Cross-correlation analysis 

Each model was convolved to the CRIRES resolution using a 
Gaussian filter, then Doppler-shifted over the range -250 to 
250kms“^ in steps of 1.5kms“\ and for each velocity step 
cross-correlated with the residual spectra. The velocity range 
was chosen to cover every possible radial velocity of the planet, 
and the step size roughly corresponds to the pixel velocity scale 
of the detectors. The outcome is a geocentric correlation matrix, 
with the strength of the cross-correlation as a function of the ap¬ 
plied radial velocity shift and the spectrum number (i.e. time). 
The cross-correlation was performed separately for each night 
and detector. 

To maximise a potential planet signal, each row of the geo¬ 
centric cross-correlation matrices, that is, the cross-correlation 
functions for each spectrum, was shifted to the rest frame of 
the planet, then summed over time. Einally, the summed, plan- 
etocentric cross-correlation functions from detectors 2 and 3 
from all three nights were combined, providing the total cross¬ 
correlation functions (T-CCE). 
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Fig. 3. Temperature-pressure profiles of models hosting a thermal in¬ 
version. The red profiles are the ones from the T/P grid described in 
Sect. |4.1| and the black profiles are the MS models described in Sect. 
|4.1.1| The deep inversion layers of the MS models give rise to strong 
CO emission. 


The radial velocity corrections (rvcon) from the geocentric to 
the planetcentric frame are a function of time and are given by 


^t’corrCO ~ ^t’pianetCO "f ^t’helioCO "f ^sys* 


( 1 ) 


Here, r^pianet is the radial velocity curve of the planet, r^heiio is the 
heliocentric correction term, and Vsys is the systemic ve locity of 


Nidever et ah 20021. Assum- 


the host star (-14.69 + 0.09 km s~ 
ing a circular orbit ( [Crossfield et al.||2012| ), the radial velocity 
curve is given by 


^i^pianet(0 = sin(27r^(0), 


( 2 ) 


with Kp being the semi-amplitude of the orbital radial velocity of 
the planet and <p being the orbital phase. The orbital period and 
the mid-transit time used to calculate the phases of the planet are 
adopted from |Knutson et al.] ( [2007b| ). 

Since HD 209458 b i s a transiting plane t with an inclination 
of i - 86.589° + 0.076° ( Southworth|2008| ) close to 90 degrees. 

equals the orbital velocity Vorb = 140 + 10kms“^ ( Snellen| 
|et al.|2010| to within the uncertainty. However, we investigated 
a wider range of values (50-170 km s“^) and Vsys values (-75 
- 45 km s“^) to scan the parameter space and check the strength 
of potential spurious signals. 


5. Results 

The diagrams in Fig. show the strength of the total cross¬ 
correlation as function of and Vsys for the four example mod¬ 
els from Fig.|^ The colour scale indicates the signal of the total 
cross-correlation function for a given parameter set. The S/N was 
obtained by dividing the total cross-correlation values by their 
standard deviation excluding points in parameter space close to 
the expected signal position. The values translate into or¬ 
bital inclinations, which are given on the right-hand axis of the 
figures. The literature and Vsys values are indicated by the 
dashed white lines. 


The absorption models and the weak-inversion models in the 
grid show a S/N = 3-3.5 hint of a CO correlation signal at a po¬ 
sition within the uncertainties of the literature values of and 
Vsys. The models from the T/P grid with the strongest emission 
lines and only low-level absorption in the wings, such as shown 
in the bottom panel of Fig. give the lowest cross-correlation 
values. In fact, these models give rise to anti-correlation sig¬ 
nals with S/N values of -2.5 —3.5. This is expected because 
the strong-inversion models are close to the inverse of the non¬ 
inversion models. The peak of the cross-correlation signal in the 
T-CCF diagram is in all cases very close (within 1 cr) to the ex¬ 
pected position of the planet signal. The values of the peak cross¬ 
correlation signals, their positions, and the cross-correlation sig¬ 
nal at the expected planet position for the four example models 
are all given in Table 

The ~ 3(T signal we obtain for non-inverted models is in our 
opinion not strong enough to claim an unambiguous detection of 
CO. In the T-CCF diagrams (Fig. 0 one can see spurious signals 
that are only slightly lower in significance. Indeed, in previous 
analyses where we claimed planet signals, we obtained signifi- 


cantly higher S/N in the range S/N 4-6 ([Snellen et al.j 

2010 

Brogi et aL|2012[|Birkby et al.|2013l|de Kok et aL|2013[ 

Brogij 

et aL|2014j). 


Table 2. Results from the cross-correlation analysis for the four exam¬ 
ple models from Fig. The first column gives the S/N of the corre¬ 
lation or anti-correlation signal at the expected planet postion (-14.69, 
140)km s"^ The second column gives the peak (anti-)correlation value, 
and the third column gives the location of the peak relative to the ex¬ 
pected position. Note that these positions are all well within the Icr 
uncertainty. 



S/N 

at planet pos. 

S/N 

peak 

peak location 
[km s“^] 

narrow abs. 

2.8 

3.2 

(-1.5, -f4.5) 

broad abs. 

3.2 

3.6 

(-f4.5, -f7.5) 

weak inv. 

3.0 

3.4 

(-f4.5, -f7.5) 

strong inv. 

-2.8 

-3.6 

(-f4.5, -f7.5) 


The diagonal noise structures in the T-CCF diagrams of Fig. 
fflare a result of a degeneracy between and Vsys, with the slope 
aepending on the observed orbital phases of the planet. The com¬ 
bination of three nights with slightly different phases, all with 
the planet on the same side of the host star, has acted to broaden 
and distort these naturally occurring noise structures. The noise 
structures are also affected by the line shape of the model. Cross¬ 
correlation with the narrow absorption model results in more nar¬ 
row noise structures, and the location of the potential CO signal 
is better defined. At the opposite end of the scale is the weak- 
inversion model, for which the line shapes at the CRIRES resolu¬ 
tion resemble absorption doublets (see Fig.[^. Cross-correlation 
with this much broader line profile smoothes the potential signal 
along the diagonal noise structures, increasing the uncertainty in 
the peak position. The peak cross-correlation values as given in 
Tableall fall along the degenerate diagonal of the diagrams. 


5.1. MS models 


The models with temperature-pressure profiles similar to the 
best-fit models of |Madhusudhan & S eager] pOQ9| ) have two basic 
types of CO line shapes at the CRIRES resolution, as described 


in Sect. 4.1.1 The pure emission lines can be seen as a mirror 
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Fig. 4. Total cross-correlation diagrams for the four models illustrated in Fig. Rafter summing over time for a range of Vsys and Kp values, and 
combining detectors 2 and 3 from the three observation nights. The dashed \^ite lines indicate the expected planet signal based on literature 
values, which coincides with a weak correlation signal in the case of the two absorption models as well as the weak-inversion model, and a weak 
anti-correlation signal in the case of the strong-inversion model. Although there are hints of a signal visible, CO is not detected with statistical 
significance. 


of the pure absorption lines from the T/P grid, and the lines with 
absorption at the core and strong emission in the wings are mir¬ 
rors of the line shapes of weak-inversion models from the T/P 
grid. We do not detect carbon monoxide when cross-correlating 
the MS models with the observed residual spectra, as described 
in Sect. |4.2| The S/N at the expected planet position is about 
-3.3 for the models with strong simple emission lines and about 
-2.8 for the models with some absorption at the core of the oth¬ 
erwise emission-dominated lines. The weak anti-correlations are 
consistent with the results from the T/P grid models, and indeed 
the T-CCF diagrams of the MS models look like colour-inverted 
versions of the diagrams for the models with broad absorption 
lines and weak inversion in Fig.|^ 


5.2. Watermodels 

We have also cross correlated model spectra with water as 
the single trace gas. The water models were calculated us¬ 
ing the same T/P grid as for the CO models, but the tested 
volume-mixing ratios for the H 2 O models were VMR = 
[10“^, 10“"^, 10“^]. We did not detect a water signal from cross- 
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correlating with any of the H 2 O models. The typical correlation 
S/N at the expected planet position lies within +1, and there is no 
systematic behaviour in the T-CCF diagrams as a fu nction of the 
line profiles, such as for the CO models. Recently, |Madhusud^ 
han et al.|(|2014|) constrained the H 2 O abundance of HD 209458 
b to sub-solar levels assuming a clear atmosphere, with the cen¬ 
tral H 2 O VMR measured to be 5.4 x 10“^, which is lower than 
the VMRs of the models in this work. The water abundance may 
still be solar or supersolar if the atmosphere of HD 209458 b has 
clouds or hazes. 


5.3. Estimating expected CO signals 

We injected model spectra into the observed spectra as fake 
planet signals and subsequently retrieved the signals to evalu¬ 
ate the expected CO signals as function of volume-mixing ratio 
and T/P profile. 

First the model planet spectrum (Fmodei) was scaled with re¬ 
spect to the stellar continuum of the parent star. We adopted the 
effect ive temperature for HD 209458 of 6092 K ( [Boyajian et al. 
2015|), and approximated the stellar spectrum from a linear inter- 
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Fig. 5. Slices through the total cross-correlation diagrams for the four atmospheric models from Fig. assuming the literature values Vsys = 
-14.69 km s“^ and - 140kms“^ The solid line is the T-CCF from cross-correlating a given model with the observed spectra, and the dashed 
line is the T-CCF from injecting the model at nominal strength into the observed spectra before telluric removal, and then cross-correlating with 
the same model. The differences between the observed and the injected T-CCFs are shown in the lower panels, and the peak signal in the difference 
T-CCFs are good estimates of expected S/N correlation signals, if the models are correct. 


polation between two Castelli-Kumcz stellar atmosphere mod¬ 
el^ with [M/H]=0.0, log g = 4.5 and the temperatures 6000 and 
6250. The scaling was subsequently performed according to 

Fsca.edO)=%^f^f, (3) 

PsW Vs I 

where Fs is the model stellar spectrum, and rp and rs are the radii 
of the planet and star. The scaled model planet spectrum was then 
convolved to match the CRIRES resolution. For each observed 
spectrum the relative radial velocity between HD 209458 b and 
the observer was calculated using Eq. (??), and the scaled model 
planet spectrum was Doppler-shifted accordingly before being 
injected into the aligned and normalised observed spectra (step 
2 in Fig.[^ using 

^obs+fake(^) — ^obs(^) ' (1 + ^scaledC^? ^)) • (4) 

We t hen proceeded with the telluric removal as described in 
Sect. |3.3| and the fake planetary signal was retrieved by cross- 
correlating with the same input model spectrum as described in 

Sect. 021 

Since the parameter Kp = 140kms“^ of the injected sig¬ 
nal is precisely known, it is sufficient to examine slices through 
the T-CCF diagram when estimating the expected signals. Such 
one dimensional T-CCfs are shown in Fig. for each of the 
four example models from Fig. The T-CCFs obtained from 
cross-correlating with the observed spectra are the solid lines, 
and the dashed lines are the T-CCFs from cross-correlating with 
the spectra injected with artificial planet signals. The injected 
planetary signatures show up as central peaks. We note that the 
shape of the four injected T-CCFs in Fig. are to a high de¬ 
gree determined by the auto-correlation function of the average 
line profile of the model. This is especially true for the relatively 

^ ftp://ftp.stsci.edu/cdbs/grid/ck04models/ 


large negative wings on either side of the correlation peak of the 
strong-inversion model. 

Since the injected signal is superimposed on top of the ob¬ 
served T-CCF, we need to examine the difference between the 
T-CCF with and without the injected signal, which are shown 
in a sub-panel below the T-CCFs in Fig. This shows that we 
can expect signals with a S/N of 2.6 for the narrow-absorption 
model, 1.9 for the broad-absorption model, 1.0 for the weak- 
inversion model, and 7.4 for the strong-inversion model. These 
results are to be interpreted as the approximate strength of the 
signal we should be able to detect if that spectral signature is 
present in the planetary atmosphere. 

Both of the absorption models have expected signals that are 
similar in strength to the observed signals, but because there are 
spurious signals of similar strength, more data are needed to con¬ 
fidently confirm the detection of CO absorption lines. For the 
strong-inversion model, it is immediately clear that the observed 
correlation results are a poor match to the expected results. The 
emission lines are expected to give rise to a strong planetary sig¬ 
nature, but there is no hint of this in the observed T-CCF. Instead, 
the observed T-CCF shows a weak anti-correlation signal, con¬ 
sistent with the possible presence of absorption lines. The weak- 
inversion model has an expected S/N < 1, revealing that it is very 
difficult to detect models with such line profiles. 

5.3.1. Line contrast 

For the simplest line profiles, those with either pure absorption 
or pure emission, it is possible to place an upper limit on the 
line contrast, that is, the line strength relative to the stellar con¬ 
tinuum. We performed this analysis for the two absorption mod¬ 
els in the top panels of Fig. The model was injected into the 
observed spectra with variable strength by multiplying Fscaied 
in Eq. ^ and Q with a variable scaling factor S. For each 
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value of S the difference between the observed and injected T- 
CCFs was determined, and we measured the line contrast of the 
model S * Fscaied(^). which gave rise to a correlation signal with 
S/N = 3. This resulted in a line contrast of ~ 3.2 x 10“ for both 
the narrow- and broad-absorption model. 


6. Discussion 


We searched for a CO signal in day side spectra of HD 209458 
b by cross-correlating with template spectra with CO as the sin¬ 
gle trace gas, and assuming many different temperature-pressure 
profiles that covered atmospheres with and without a tempera¬ 
ture inversion. The presence of an inversion layer gives rise to 
CO emission in the spectra of the planet, but depending on the 
pressure range and temperature gradient of the inverted layer, 
the shape of the CO lines can differ significantly. A thermal in¬ 
version will give rise to strong emission when it is located deep 
in the atmosphere and has a large temperature difference. The 
depth of a temperature inversion depends on the opacity sources 
responsible for the inversion, but since the sources are unknown 
in hot Jupiters, it is relevant to test models with inversions at dif¬ 
ferent depths. At the high spectral resolution of the CRIRES in¬ 
strument, more complex line shapes can be partly resolved, with 
the core of the lines probing higher altitudes, while the wings 
of the lines probe deeper, just above where the atmosphere be¬ 
comes optically thick. The line shape for a given T/P profile is 
to a lesser degree also affected by the CO volume-mixing ratio 
because the VMR influences the probed pressure range. 

We have tested models with a wide variety of line shapes: 
Pure absorption and pure emission lines, as well as models with 
either emission at the core and absorption in the wings, or ab¬ 
sorption at the core and emission in the wings. We did not de¬ 
tect carbon monoxide in the HD 209458 b day side spectrum 
with statistical significance, regardless of which model we cross- 
correlated with. The models with non-inverted atmospheres and 
the models with a weak-inversion layer, typically located at pres¬ 
sures lower than 0.01 bar, all show a hint of a correlation signal 
with S/N = 3.0 - 3.5, while the models with stronger inver¬ 
sion layers, including the MS models with T/P profiles like those 
of |Madhusudhan & Seager|p009| ), show anti-correlation signals 
with similar low-significance levels. 

Our results are inconsistent with the inference of IKnutsonI 
|et al.| ( [200^ of a thermal inversion in HD 209458 b based on 
Spitzer IRAC broadband photometry, but agrees with the more 
recent IRAC data and anal ysis by [Zellem et al.| ( |2014| ) and 
|Diamond-Lowe et ar| ( |2014| ), who both found that a thermal in¬ 
version is not required to explain the Spitzer secondary-eclipse 
depths. In fact, the best-fitting T/P profile found by [Diamond^ 
|Lowe et ak] ( |2014| ) is isothermal at pressures lower than 0.1 
bar. However, the high-resolution spectroscopy method is only 
sensitive to narrow components in the planetary spectrum - any 
broad component is filtered out in the early stages of the analy¬ 
sis. In Fig.we show the contribution function for a CO volume¬ 
mixing ratio of 10“^, indicating that our data are mostly sensitive 
to temperature differences in the 10“^ to 10“^ pressure range. 
This means that we expect no signal from an atmosphere such 
as shown by Diamond-Lowe et al.|p014] ). Our non-detection is 
therefore consistent with their results. 

The non-detection of CO in the dayside spectrum of HD 
209458 b is interesting in light of previous successes with de¬ 
tecting CO using high-resolution spectrosco py at 2.3 pm, both 
in the dayside spectra of other hot J upiters ([Brogi et al^2012[ 


Rodler et al.|2012[|Brogi et al.|2013t[de Kok et al.|2013[ Brogi| 


et al.||2014| ) and in the transmission spectrum of HD 209458 b 
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Fig. 6. Contribution function of the CRIRES observations for a CO 
volume-mixing ratio of 10"^, with CO as a single trace gas. Note that 
these observations are sensitive to the narrow-line components over a 
pressure range -10~^ to -10~^ bar. In particular, i n this range the T/P 
profile as derived by |Diamond-Lowe et"^ ( |2014| ) is near isothermal, 
making these observations insensitive to it. 
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itself ( [Snellen et al.|2010| ). [Snellen et al. [ detected CO absorption 
and determined the line strength relative to the stellar continuum 
to be 1 - 1.5 X 10“^, while in this work we obtain a best esti¬ 
mate of the relative line strength for simple absorption lines in 
the dayside spectrum to be ~ 3.2 x 10“"^. 

The most likely explanation for the lack of CO signal is that 
the average dayside T/P profile is near-isothermal in the pressure 
range we probe (10“^ - 10“^ bar, see Fig.[^. However, we note 
that at high altitudes there may still be a weak inversion layer, 
for instance at pressures « 10“^ bar, which are only probed by 
very narrow cores of the CO lines. Furthermore, when observ¬ 
ing the dayside spectrum of HD 209458 b, the spectrum is disk 
integrated. HD 209458 b is expected to be tidally locked to its 
parent star, and [Zellem et al.[ ( [2014[ ) have observed the planet 
to have a hotspot, shifted eastward of the substellar point due 
to equatorial superrotation, as expected if the planet is in syn¬ 
chronous rotation. It is possible that the planetary atmosphere 
has a thermal inversion only in the substellar point, or alterna¬ 
tively, in the eastward-shifted hotspot. If a thermal inversion is 
localised to only a part of the dayside, the effect of the inversion 
will be muted. 

Although high, the limited spectral resolution (R = 100000) 
of the CRIRES instrument can work to dampen a thermal in¬ 
version signal somewhat. This is especially true for the weak¬ 
est inversions (high altitude and less steep temperature gradient 
dT/dP), where the emission at the core of the lines is so weak 
that the complex line profiles are reduced to slightly muted, sim¬ 
ple absorption lines. This effect alone cannot explain the lack of 
a CO detection, but it may be a contributing factor. 

Alternatively, the CO signal, whether it is absorption or emis¬ 
sion, might be heavily dampened by clouds or hazes. [Charbon- 
[neau et ( [2002[ ) first suggested an opaque cloud layer as a 
possible explanation for the lack of a strong optical Na feature 
in the optical transmission spectrum of HD 209458 b, although 
they noted that if clouds were the sole explanation, the cloud 
tops were required to be located above pressures of 0.4 mbar. 
It is unclear if such high-altitude cloud layers can exist in hot 
Jupiters, although [Fortney et al.[ ( [2003) suggests it is possible to 
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have clouds at pressures of several millibars. However, one does 
expect that the contributions from clouds or hazes will be signif¬ 
icantly stronger in transit spectra than in the dayside spectrum 
becaus e of the grazing geometry in th e former. [Deming et al.| 
( |2013| ) and |Madhusudhan et al.| ( |2Ql^ have observed the trans¬ 
mission spectrum of HD 209458 b with HST/WFC3 and both 
found weak transit features that possibly require haze. The haze 
alone is not expected to have a major influence on the thermal 
emission spectrum, but a cloud or haze layer at altitudes below 
those probed by transit observations might influence the contin¬ 
uum level for the thermal emission spectra and hence our results. 
The role of clouds or haze is thus still unclear. 


7. Conclusion 

We did not detect carbon monoxide at a statistically signifi¬ 
cant level in the 2.3 qm high-resolution dayside spectrum of HD 
209458 b, although we saw a potential absorption signal at the 3 
- 3.5(T level. CO is e xpected to be abunda nt in hot Jupiters ac¬ 
cording to theory (e.g. |Madhusudhan|2Q12| ), and has been found 
to be present in the tranmission spectrum of the planet ( |Snellen| 
|et al.|2010| ). A near-isothermal pressure profile of the planetary 
atmosphere is the most probable cause of the lack of CO signal. 

We showed that high-resolution spectroscopy can in princi¬ 
ple measure line shapes and constrain the temperature pressure 
profiles in exoplanet atmospheres. In general, it will be possible 
to distinguish between a clear atmosphere with and without a 
thermal inversion, as long as these occur within the probed pres¬ 
sure range. 

Our observations support the recent findings by |Zellem et al. 
( 2014| ) and |Diamond-Lowe et al.| ( [2014| ) that there is no strong 
inversion layer in the dayside atmosphere of HD209458b, in the 
pressure range 10“^ bar - 10“^ bar. The emission signal from 
such an inversion layer would have been readily detected in our 
spectra. 

Acknowledgements. We are thankful to the ESO staff of Paranal Observatory for 
their support during the observations, and we thank the anonymous referee for 
the insightful comments and suggestions. This work is part of the research pro¬ 
grammes PEPSci and VICI 639.043.107, which are financed by the Netherlands 
Organisation for Scientific Research (NWO). Support for this work was provided 
in part by NASA, through Hubble Fellowship grant HST-HF2-51336 awarded by 
the Space Telescope Science Institute. This work was performed in part under 
contract with the California Institute of Technology (Caltech)/Propulsion Labo¬ 
ratory (JPL) funded by NASA through the Sagan Fellowship Program executed 
by the NASA Exoplanet Science Institute. 


References 

Arsenault, R., Alonso, J., Bonnet, H., et al. 2003, in Adaptive Optical System 
Technologies II, Vol. 4839, SPIE, 174 
Barnes, J. R., Barman, T. S., Prato, L., et al. 2007, MNRAS, 382, 473 
Bento, J., Wheatley, P. J., Copperwheat, C. M., et al. 2014, MNRAS, 437, 1511 
Birkby, J. L., de Kok, R. J., Brogi, M., et al. 2013, MNRAS, 436, L35 
Borysow, A. 2002, A&A, 390, 779 

Borysow, A., Jprgensen, U. G., & Fu, Y. 2001, J. Quant. Spec. Radiat. Transf., 
68, 235 

Boyajian, T., von Braun, K., Feiden, G. A., et al. 2015, MNRAS, 447, 846 
Brogi, M., de Kok, R. J., Birkby, J. L., Schwarz, H., & Snellen, I. A. G. 2014, 
A&A, 565, A124 

Brogi, M., Snellen, I. A. G., de Kok, R. J., et al. 2012, Nature, 486, 502 
Brogi, M., Snellen, I. A. G., de Kok, R. J., et al. 2013, ApJ, 767, 27 
Brown, T. M., Charbonneau, D., Gilliland, R. L., Noyes, R. W., & Burrows, A. 
2001, ApJ, 552, 699 

Brown, T. M., Libbrecht, K. G., & Charbonneau, D. 2002, PASP, 114, 826 
Burrows, A., Budaj, J., & Hubeny, I. 2008, ApJ, 678, 1436 
Burrows, A., Hubeny, L, Budaj, J., Knutson, H. A., & Charbonneau, D. 2007, 
ApJ, 668, L171 

Charbonneau, D., Allen, L. E., Megeath, S. T., et al. 2005, ApJ, 626, 523 


Charbonneau, D., Brown, T. M., Latham, D. W., & Mayor, M. 2000, ApJ, 529, 
L45 

Charbonneau, D., Brown, T. M., Noyes, R. W., & Gilliland, R. L. 2002, ApJ, 
568,377 

Charbonneau, D., Knutson, H. A., Barman, T., et al. 2008, ApJ, 686, 1341 
Christiansen, J. L., Ballard, S., Carbonneau, D., et al. 2010, ApJ, 710, 97 
Crossfield, I. J. M., Knutson, H., Fortney, J., et al. 2012, ApJ, 752, 81 
de Kok, R. J., Brogi, M., Snellen, I. A. G., et al. 2013, A&A, 554, 82 
de Pater, I. & Lissauer, J. J. 2010, Planetary Sciences, 2nd edn. (Cambridge 
University Press), ISBN: 978-0-521-85371-2 
Deming, D., Brown, T. M., Charbonneau, D., Harrington, J., & Richardson, L. J. 
2005, ApJ, 622, 1149 

Deming, D., Seager, S., Richardson, L. J., & Harrington, J. 2005, Nature, 434, 
704 

Deming, D., Wilkins, A., McCullough, P, et al. 2013, ApJ, 774, 95 
Desert, J.-M., Vidal-Madiar, A., Lecavelier Des Etangs, A., et al. 2008, A&A, 
492, 585 

Diamond-Lowe, H., Stevenson, K. B., Bean, J. L., Line, M. R., & Fortney, J. J. 
2014, ArXiv 1409.5336 

Fortney, J. J., Fodders, K., Marley, M. S., & Freedman, R. S. 2008, ApJ, 678, 
1419 

Fortney, J. J., Saumon, D., Marley, M. S., Fodders, K., & Freedman, R. S. 2006, 
ApJ, 642, 495 

Fortney, J. J., Sudarsky, D., Hubeny, I., et al. 2003, ApJ, 589, 615 
Gibson, N. P, Aigrain, S., Barstow, J. K., et al. 2013, MNRAS, 436, 2974 
Guillot, T. 2010, A&A, 520, All 

Hoeijmakers, H. J., de Kok, R. J., Snellen, I. A. G., et al. 2014, ArXiv e-prints 
Horne, K. 1986, PASP, 98, 609 

Hubeny, I., Burrows, A., & Sudarsky, D. 2003, ApJ, 594, 1011 
Hudson, C. M., Sing, D. K., Pont, F, et al. 2013, MNRAS, 434, 3252 
Kaufi, H., Ballester, P, Biereichel, P, et al. 2004, in Ground-based Instrumenta¬ 
tion for Astronomy, Vol. 5492, SPIE, 1218 
Knutson, H. A., Charbonneau, D., & Allen, L. E. 2008, ApJ, 673, 526 
Knutson, H. A., Charbonneau, D., Allen, L. E., et al. 2007a, Nature, 447, 183 
Knutson, H. A., Charbonneau, D., Burrows, A., O’Donovan, F. T., & Mandushev, 
G. 2009, ApJ, 691, 866 

Knutson, H. A., Charbonneau, D., & Noyes, R. W. 2007b, ApJ, 655, 564 
Knutson, H. A., Howard, A. W., & Isaacson, H. 2010, ApJ, 720, 1569 
Line, M. R., Knutson, H., Wolf, A. S., & Yung, Y. L. 2014, ApJ, 783, 70 
Lockwood, A. C., Johnson, J. A., Bender, C. F., et al. 2014, ApJ, 783, L29 
Lord, S. D. 1992, A New Software Tool for Computing Earth’s Atmospheric 
Transmission of Near- and Far-infrared Radiation, Technical Memorandum 
103957, NASA 

Machalek, R, McCullough, P. R., Burke, C. J., et al. 2008, ApJ, 684, 1427 
Machalek, R, McCullough, P. R., Burrows, A., et al. 2009, ApJ, 701, 514 
Madhusudhan, N. 2012, ApJ, 758, 36 

Madhusudhan, N., Crouzet, N., McCullough, P. R., Deming, D., & Hedges, C. 
2014, ApJ, 791, L9 

Madhusudhan, N., Mousis, O., Johnson, T. V., & Lunine, J. 1. 2011, ApJ, 743, 
191 

Madhusudhan, N. & Seager, S. 2009, ApJ, 707, 24 
Madhusudhan, N. & Seager, S. 2010, ApJ, 725, 261 

Moses, J. I., Madhusudhan, N., Visscher, C., & Freedman, R. S. 2013, ApJ, 763, 
25 

Nidever, D. L., Marcy, G. W., Butler, R. P, Fischer, D. A., & Vogt, S. S. 2002, 
ApJS, 141,503 

Oberg, K. I., Murray-Clay, R., & Bergin, E. A. 2011, ApJ, 743, L16 
O’Donovan, F. T., Charbonneau, D., Harrington, J., et al. 2010, ApJ, 710, 1551 
Parmentier, V, Showman, A. R, & Lian, Y. 2013, A&A, A91 
Rodler, F., Lopez-Morales, M., & Ribas, I. 2012, ApJ, 753, L25 
Rothman, L. S., Gordon, I. E., Barber, R. J., et al. 2010, J. Quant. Spec. Radiat. 
Transf., 111,2139 

Schneider, J. 1994, Ap&SS, 212, 321 

Seager, S. 2010, Exoplanet Atmospheres: Physical Processes (Princeton Univer¬ 
sity Press), ISBN: 978-0-691-13027-9 
Seager, S. & Deming, D. 2010, ARA&A, 48, 631 
Seager, S. & Sasselov, D. D. 2000, ApJ, 537, 916 

Sing, D. K., Lecavelier des Etangs, A., Fortney, J. J., et al. 2013, MNRAS, 436, 
2956 

Snellen, I. A. G., de Kok, R. J., de Mooij, E. J. W., & Albrecht, S. 2010, Nature, 
465,1049 

Snellen, I. A. G., de Kok, R. J., le Poole, R., Brogi, M., & Birkby, J. 2013, ApJ, 
764,182 

Snellen, I. A. G., de Mooij, E. J. W., & Albrecht, S. 2009, Nature, 459, 543 
Southworth, J. 2008, MNRAS, 386, 1644 
Spiegel, D. S., Silverio, K., & Burrows, A. 2009, ApJ, 699, 1487 
Webb, J. K. & Wormleaton, I. 2001, PAS A, 18, 252 

Zahnle, K., Marley, M. S., Freedman, R. S., Fodders, K., & Fortney, J. J. 2009, 
ApJ, 701, L20 

Zellem, R. T., Lewis, N. K., Knutson, H. A., et al. 2014, ApJ, 790, 53 


Article number, page 11 ofpT] 



























